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Abstract. Optical reference geometry and related concept of inertial forces are 
investigated in Kerr-de Sitter spacetimes. Properties of the inertial forces are 
(-H ' summarized and their typical behaviour is illustrated. The intuitive 'Newtonian' 

' application of the forces in the relativistic dynamics is demonstrated in the case 

of the test particle circular motion, static equilibrium positions and perfect fluid 
toroidal configurations. Features of the optical geometry are illustrated by the 
embedding diagrams of its equatorial plane. The embedding diagrams do not 
cover whole the stationary regions of the spacetimes, therefore the limits of 
embeddability are established. A shape of the embedding diagrams is related 
^ ' to the behaviour of the centrifugal force and it is characterized by the number 

of turning points of the diagrams. Discussion of the number of embeddable 
s ! ■ photon circular orbits is also included and the typical embedding diagrams 

C D i are constructed. The Kerr-de Sitter spacetimes are classified according to the 

properties of the inertial forces and embedding diagrams. 
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1. Introduction 



In 1988, a new approach to the relativistic dynamics appeared. By an appropriate 
conformal (3+1) splitting of the Schwarzschild spacetime, the so-called optical 
reference geometry was defined pQ. It was demonstrated that this geometry enables 
us to define the inertial forces in the framework of general relativity, providing a 
description of test particle motion in accordance with our 'Newtonian' intuition [2] , [3] . 
Later, the definition of the optical reference geometry and related inertial forces was 
generalized, extending its applicability to any spacetime [I], [5], [5J, [7J. The optical 
reference geometry was also thoroughly studied in particular spacetimes, such as 
Schwarzschild-de Sitter [B] , [9] , Reissner-Nordstrom [TO] , Reissner-Nordstrom-de Sitter 
[11], Kerr [12], [13] or Kerr- Newman [14], illustrated by embedding diagrams [15] . 
|16] , [17j , and the inertial forces formalism was applied for solving specific problems in 
relativistic dynamics [18] , [19] . Behaviour of the centrifugal force is closely related to 
the shape of embedding diagrams of the optical geometry, therefore many properties 
of the relativistic dynamics in the spacetimes were effectively illustrated and directly 
visualized. 

In the present work we extend our previous studies [9], [13] to the case of the 
stationary and axially symmetric Kerr-de Sitter (KdS) spacetimes, reflecting thus 
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their basic properties incorporating a combined influence of the rotation of source 
and the cosmic repulsion, recently indicated by wide range of cosmological tests. 
It is worth to stress that the KdS black-hole spacetimes could be important in 
understanding astrophysical phenomena exposed around supermassive black holes 
in giant active galactic nuclei as demonstrated in [20], [21]. We show that the 
application of the inertial forces concept proves to be very effective and provides 
simpler and intuitive approach to some relativistic problems, as compared to the 
standard general relativistic methods. The standard methods based on the covariant 
formulation of physical laws in a given background (geodesic equation, energy 
momentum conservation, etc.) can be then replaced by the well developed methods 
familiar from Newtonian physics when inertial forces related to the optical geometry 
are used, making a significant technical simplification in treating physical processes 
in curved backgrounds and enabling their interpretation in the Newtonian way. Thus 
the present paper could be also considered to be an introduction in developing a new 
approach to study astrophysical phenomena. We give some explicit examples how the 
analysis of the inertial forces could immediately and in a very simple way make an 
enlightenment of the properties of the circular motion of test particles and perfect 
fluid. 

The structure of the paper is similar to the previous study of the optical geometry 
of the Kerr-Newman spacetimes |14j . whereas the content is more extensive, mainly 
in the parts concerning the inertial forces application in the relativistic dynamics. 
In section 2, we briefly present basic general definitions of the optical geometry and 
inertial forces, focusing on the stationary and axially symmetric spacetimes and 
circular motion. In section 3, KdS geometry is introduced. Although the properties 
of the KdS spacetimes are well known, we briefly discuss them here again, in order 
to simply demonstrate the so-called 'Chinese boxes' method, used later in much more 
complicated cases for discussions of the properties of the inertial forces and embedding 
diagrams. In section 4, the inertial forces in the KdS spacetimes are determined for 
general circular motion, and their properties are briefly discussed and thoroughly 
represented for the equatorial circular motion. In section 5, some examples of the 
relativistic dynamics are shown from the point of view of the inertial forces formalism. 
There are many interesting relativistic problems, which can be simply investigated by 
using the inertial forces formalism. We present such problems that have been studied 
by the standard general relativistic methods as well, in order to show immediately 
the effectiveness and simplicity of the inertial forces approach. Namely the equatorial 
circular motion is studied, simply providing important findings for the test particle and 
photon circular motion [22], equatorial static radius and static equilibrium positions 
of test particles on the axis of symmetry [23 . Another astrophysically relevant 
application of the inertial forces formalism is demonstrated in deducing projection 
of the relativistic Euler equation that allows determination of the perfect fluid 
equilibrium toroidal configurations [20] . In section 6, the embedding diagrams of the 
equatorial plane of the optical reference geometry of KdS spacetimes are discussed and 
constructed. The whole optical geometry cannot be embeddable, therefore we discuss 
limits of embeddability. We establish relations between the shape of the embedding 
diagrams and the behaviour of the centrifugal force, investigating turning points of 
the diagrams. Finally, the embeddability of photon circular orbits is discussed and 
typical embedding diagrams are constructed. In section 7, we summarize results of 
our investigation of the optical reference geometry formalism in the KdS spacetimes, 
and compare them with results obtained for previously studied backgrounds. 
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2. Fundamental definitions of the formalism 

2.1. General case 

In stationary spacetimes described by a metric gik (with signature +2), the definition 
of the optical reference geometry requires an introduction of family of special observers 
with a timelike, unit, and hypersurface orthogonal 4-velocity field n l and with its 
4-acceleration field equal to the gradient of a scalar function, i.e., [5] 

n k n k = -1, n l V l n k = V fc $, n^VjU).] = 0. (1) 

Such a vector field can be chosen in the form 

rj=e-*L\ $ = iln(-t < 4 i ), (2) 

whereas it corresponds to the unit 4-velocity field of stationary observers parallel to 
a timelike vector field l 1 . Then the local instantaneous 3-dimensional space of the 
observers is described by the metric 

= gik + riiUk, (3) 

the so-called directly projected geometry. Its conformally adjusted metric 

kk = e~ 2 *h lk (4) 

is the so-called optical reference geometry. 

The 4-velocity u l of a test particle with a rest mass m can be decomposed in the 
reference frame of the special observers with the 4-velocity n % by using the relation 

u* =-y(n i + VT i ), (5) 

where r l is the unit spacelike vector parallel to the 3-velocity v % of the particle in the 
3-dimensional space and 7 = (1 — v 2 )~ x / 2 is the Lorentz factor (the speed of light 
c = 1). Moreover, there is 7 = —tfui, jvt 1 = 71/ = u k h k , where h\ = 5\ + n l nk is a 
projection tensor. Note that the vector t % is defined only along the world line of the 
particle. But in our construction, it is necessary to know how r 4 changes along n l . 
This can be done in several different ways (different gauges). The gauge used by M. 
Abramowicz and his co-workers assures that the Lie derivation of t 1 with respect to 
L l = e*n l vanishes, i.e., C L Tk = t'VjTfc — r'Vitfc = 0, which causes that the Coriolis 
force defined below vanishes in static spacetimes [5]. 

The 4-acceleration ak = u l ViUk of the particle can be written in the form 

a k = 7 2 V fc $ + 7 2 v(n l V 2 T fc + T l V t n k ) + j 2 v 2 T l V z T k + {vj) Tk + 7K/6) 

where (vj) = u l V i{^fv). By using the spacelike unit vector parallel to t % in the optical 
reference geometry , i.e., the vector f* = e*r l , its covariant form fj = e~*Ti, the 
scalar E — —i l v,i, identity j 2 = 1 + v 2 j 2 and condition of hypersurface orthogonality, 
the projection of the 4-acceleration into the 3-dimensional space can be written as 

af = h)ak - Vi$ + ( 7 t0 2 f l V 4 f, + l 2 vX, + Vfj, (7) 

where Xj = n 4 (ViT,- — Vjn) and V — u' L Vi(Ev). 

Now, the inertial force in the 3-dimensional space -F^" 1 ", related to the real force 
by the equation Fj, 1 - = — F^, can be decomposed into the sum of the gravitational G k , 
centrifugal Zk, Coriolis CV- and Euler Ek forces familiar from the Newtonian physics, 
i.e., 

Fk ± = -mai = G k + Z k +C k + E k , (8) 
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where 

G k = -mVk$, (9) 

Z k = -m( 7U ) 2 f l V l f fe , (10) 

C k = ~ mi 2 vX kl (11) 

E k = -mVf k . (12) 

Note that in the context of the optical reference geometry approach, the gravitational 
force ranks among the inertial forces. 

2.2. Axially symmetric spacetimes and circular motion 

In static and spherically symmetric spacetimes, the 4-velocity held ([2]) can be chosen 
with the vector field i % being the timelike Killing vector field rf = 5\, which is 
orthogonal to the spacelike azimuthal Killing vector field £ l = 8\. However, in 
stationary and axially symmetric spacetimes, the Killing vector fields rf and £ l are 
not orthogonal in general. Thus the vector field i % has to be chosen in order to satisfy 
the hypersurface orthogonality condition. It is easy to check that the vector field 
l 1 = rf + flLNRF^j where Qlnrf = —rfii/di IS hypersurface orthogonal and thus 
the 4-velocity field 

n l ^e-^tf + n LNRF e), (13) 

$ = ^ In [-(rf + VL LNRF C){m + SILNHF&)], (14) 

which corresponds to the 4-velocity field of the locally non-rotating frames moving 
along circular orbits with the angular velocity d<j>/dt = Qlnrf, defines the special 
observers [5]. Since the vector field n 1 is hypersurface orthogonal, the unit vector r 4 , 
used in the decomposition ([5]), must be located in the hypersurface. On the other 
hand, the 4-velocity 

u* = e-'V + ( 15 ) 

A=Un[-(r 1 i + ne)(Vi + ^i)], (16) 

corresponding to the general circular motion of a test particle with an angular velocity 
d<j)/dt = Q =/= Qlnrf, is not hypersurface orthogonal. Moreover, the circular motion 
is directed along the Killing vector i.e., 

r i = (t%r 1/2 e- (17) 

Then, by using expressions ([5]) and (|15p . the Lorentz factor and velocity are given by 
the relations 

7 = e*~ A , v = e-*(Z%) 1 / 2 (n-n LNRF ) (18) 
and components of the inertial forces (f9j)- ([T2"l) are given by the relations 

G fe = -mVfc$, (19) 

Z k = m{ivf i(f6)^e- 2 *[e 2 *V fc (f &) - f &V fe e 2 *], (20) 

C k = mi 2 v ^%r 3/2 e-^[C^k(v^i) - V%Vk(?&], (21) 
E k =me^ s u i V i {v)f k . (22) 
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The 'Newtonian' character of the inertial forces formalism appears after introducing 
the quantities 

£=(e&) 1/2 e-*, n = n-n LNRF , v = lv , (23) 

where f2 is angular velocity relative to the LNRF. Then we can write v = Rti and the 
inertial forces IT9 |) - ([22|) can be rewritten to the form 

G k = -mV fc $, (24) 

Z k = mv 2 R- l V k R, (25) 

C k = -m(l + v 2 ) 1 / 2 vRV k ^LNRF, (26) 

E k = m(l + i 2 f /2 e*fl« < V ( (n)^. (27) 

The quantity i? corresponds to the so-called radius of gyration, defined by the relation 

R = JJL (28) 
V Eil y ' 

generalizing in a special way its definition introduced and discussed for the static 
Schwarzschild spacetimes in [26], [5]. The quantity L — u l £i is the conserved 
specific angular momentum and E = —u l ii is a scalar called 'pseudo-energy' of 
the particle. Since the vector field l 1 is not a Killing vector field, such 'pseudo- 
energy' is not conserved in general stationary and axially symmetric spacetimes. Of 
course, in the static spacetimes, E is conserved, being defined in the standard way 
as E = E = —u l r\i. The radius of gyration plays an important role in the theory of 
rotational effects in strong gravitational fields. The direction of increasing radius of 
gyration gives a preferred determination of the local outward direction relevant for 
the dynamical effects of rotation, whereas the direction becomes misaligned with the 
'global' outward direction in strong fields [26] . 

3. Kerr-de Sitter spacetimes 

In the standard Boyer-Lindquist coordinates x l — (t, r, 9, <p) and geometric units 
(c=G=l), the line element of the KdS geometry is given by the relation 

ds 2 = -£f(dt -asm 2 Odcji) 2 + Ae T f^ 9 [adt - (r 2 + a 2 )#] 2 + 

+ (29) 



where 



A r = r 2 ~2Mr + a 2 ~--Ar 2 {r 2 + a 2 ), (30) 



3 

Ag = 1 + \ka 2 cos 2 9, (31) 



1=1 + ^Aa 2 , (32) 

p 2 = r 2 + a 2 cos 2 6. (33) 

The mass M, specific angular momentum a, and cosmological constant A are 
parameters of the spacetime. It is convenient to introduce a dimensionless cosmological 
parameter 

y = ^AM 2 . (34) 
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For simplicity, we put M = 1 hereafter. Equivalently, also the coordinates t, r, 
the line element ds, and the parameter a are expressed in units of M and become 
dimensionless. 

The KdS spacetimes, being stationary and axially symmetric, admit both the 
Killing vector fields rf and £\ whereas rfr]i = g u , = g<f><f> and rf^i = g t< p. The 
only intrinsic singularity of the KdS solution is the ring singularity in the equatorial 
plane and it is given by the relation p — 0. Stationary regions of the spacetimes are 
determined by the relation A r (r; a 2 , y) > and limited by the inner and outer black- 
hole horizons at r^~ and Th+ and by the cosmological horizon at r c , which are real roots 
of the equation A r (r; a 2 ,y) — 0. Then spacetimes containing three horizons are black- 
hole spacetimes, while spacetimes containing one horizon (the cosmological horizon 
exists for any choice of the spacetime parameters) arc naked-singularity spacetimes. 
Spacetimes with two horizons are extreme black-hole or extreme naked-singularity 
spacetimes [22] . 

In order to obtain the number of horizons in dependence on the rotational and 
cosmological parameters a 2 and y, we use the 'Chinese boxes' technique that will 
be used later in more complex discussions of the properties of the inertial forces 
and embedding diagrams. In the KdS spacetimes, functions of the radius and two 
spacetime parameters are usually considered, and the routine is realized in two 
steps. A property of a function can be examined by investigation of properties of 
another function of the radius and of the lower number of parameters, i.e., of the 
only spacetime parameter. These properties of the function can be examined by 
investigation of properties of other functions, but now, functions of the radius only. 
Thus the properties of these functions can be determined by using a common way of 
investigation of functions dependent on one variable. Finally, the KdS spacetimes can 
be classified according to the initial property of the initial function. 

In this case, it follows from the relation A r (r;a 2 ,y) = that for given a 2 and 
y > (repulsive cosmological constant), the loci of horizons are given by solutions of 
the equation 

a _ r 2 - 2r + a 2 
y = y h {r; a ) = r2(r , 2 + Q2) ■ (35) 

The reality condition A r (r; a 2 ,y) > can be then rewritten in the form 

0<y<y h (r;a 2 ). (36) 

The asymptotic behaviour of the function yh(r; a 2 ) is given by y/t(r — > oo; a 2 ) — > +0 
and yh(r — * 0;a 2 ) — > oo. The local extrema of yh(r;a 2 ) are determined (due to the 
condition d r yh(r; a 2 ) = 0) by the relation 

a 2 (r)=a 2 he (r) = ^(-2r 2 + V8r~+lr + r), (37) 

the maximum of the function a^ e (r) is located at r = 1.6160 and takes the value 
a he max = 1-2120 (see figureQJ,). We can distinguish three different types of behaviour 
of yh{r; a 2 ) (see figure [TJ:). 

• For a 2 < a 2 hemax , y h (r\a 2 ) has two local extrema y h , min (a 2 ) and y h , ma x(a 2 ) 
determined by relations (|33|) and ([57|) (see figure [Tb). The black- hole spacetimes 
exist for yh,min{<i 2 ) < y < yh,max(a 2 ) and y > 0, while the naked-singularity 
spacetimes exist for < y < yh,min{a 2 ) or y > yh,rnax(a 2 ). 

• For a 2 = al e max , the extrema yh,min(a 2 ) and y h ,max(a 2 ) coincide at y h , cr it = 
0.0592, which is the limiting value for the black-hole spacetimes. 
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Figure 1. (a) Characteristic function a\ {r) governing the location of extrema 
of the function yh(r; a, 2 ). For given a 2 , the extrema of j/^(r;a 2 ) are determined 
by solutions of a 2 = a? (r) (note the dashed line a 2 = 1.06 and compare with 
figure [Tp). (b) Function ^(rja 2 ) determining the loci of event horizons of the 
KdS spacetimes and limiting the dynamic region where y > j/h(r;a 2 ) (gray). 
The function is given for a 2 = 1.06. For given a and y > 0, the horizons 
are determined by solutions of y = ?/;j(r;a 2 ). In the parameter line (y) of the 
spacetimes, the extrema of yh{r\ a 2 ) separate regions corresponding to the black- 
hole (BH) spacetimes and naked-singularity (NS) spacetimes for given values 
of the rotational parameter a 2 , (c) Functions yh(r;a 2 ) in three specific cases; 
a 2 = 0.9 (with one positive and one negative extrema), a 2 = 1.05 (with two 
positive extrema) and a 2 = 1.4 (with no extrema). (d) Functions y h m i„(a 2 ) and 
yh,max( a2 ) separating the parametric plane (a 2 ,y) into regions corresponding to 
the BH and NS spacetimes. 



• For a 2 > cbf lemax i yh{f",a 2 ) has no extrema and there are only the naked- 
singularity spacetimes for any positive value of y. 

The parameter plane (a 2 , y) separated by the functions yh,min(ci 2 ) and yh,max{o?) 
into the regions corresponding to the black-hole and naked-singularity spacetimes is 
illustrated in figure QJ1. 

4. Behaviour of inertial forces for circular motion 

Assuming a test particle circular motion with the 4-velocity (JTSJ) at a fixed radius 
and latitude in the reference frame x l = (t, r, 9) in the stationary and axially 
symmetric KdS spacetimes, there are always zero time and azimuthal components 
of the gravitational, centrifugal and Coriolis forces (|19p - (pil"]) . and zero radial and 
latitudinal components of the Euler force ([22]) . The centrifugal, Coriolis and Euler 
forces vanish for v = and moreover, the Euler force also vanishes in the case of the 
uniform motion v = const. However, due to the behaviour of the velocity independent 
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parts of the forces, they can vanish (and change their sign) at some radii and latitudes 
independently of the velocity. The same can happen for the gravitational force which 
is velocity independent by definition. Note that we pay no more attention to the 
Euler force, because of its vanishing in the case of the uniform motion considered in 
the following. Then it is convenient to rewrite the inertial forces to the form 

G k =mg k , Z k = mi, 2 Z k , C k = mv(l + v 2 ) 1 ' 2 C k , (38) 

whereas the non-vanishing velocity independent parts of the accelerations are given 
by the relations 



rA r + A A g (AB + 2rA r )p 2 
VA r p 2 



Gr = > ( 39 ) 



_ a 2 sin 2 9A r (2rA r + Bp 2 ) + AA e [2r{2p 2 - A)A r + ABp 2 } 
r ~ VA^ 2 ' ( 0) 

- _ 2asm9[rA 2 A g - ABp 2 - rA r (A + p 2 )}VA^ 

Lr ~ VA rP 2 ' [ ' 

a 2 sm9cos9(VAg-IA r p 2 ) 
= ^r— 5 , (42) 



VAgp 2 

cot 8[2AVAg - I(a 2 A e + a 2 sin 2 9A r )p 2 s 
VAgp 2 

4a 3 r cos 9 sin Q^J A r Ag 

VAgp 2 



where 



A = r 2 + a 2 , B=l-r + yr(a 2 + 2r 2 ), (45) 
V = A 2 A g - a 2 A r sin 2 6. (46) 

The presented components of the accelerations are considered to be functions of two 
variables r and 6, and two spacetime parameters a and y. There is T> > and 
Ag > for all the values of r and 9. Thus the functions are well defined throughout 
the stationary regions. The functions diverge at the radii of event horizons and in 
the ring singularity in the equatorial plane. Moreover, the latitudinal component of 
the velocity independent part of the centrifugal acceleration diverges at the axis of 
symmetry (9 — 0). 

In the case of the equatorial motion (8 = ir/2), latitudinal components of all 
the forces vanish. For given values of the rotational and cosmological parameters 
a 2 and y > 0, the radii of circular orbits where the velocity independent parts of 
the radial components of the gravitational, centrifugal and Coriolis accelerations, 
i.e., G r (r;a 2 ,y), Z r (r; a 2 ,y) and C r (r;a 2 ,y), vanish are given by zero points of the 
numerators of relations (|3"5j) - (|4*Tj) , provided that y < yh(r;a 2 ) because of the reality 
condition. It is immediately clear from relation (|4Tj) that C r (r; a 2 , y) does not vanish at 
all. Moreover, since the zero points of the functions Q r (r\ a 2 ,y) and Z r (r; a 2 ,y) enable 
changes of their sign, the gravitational and centrifugal forces also can change their 
orientation on such orbits. Detailed discussion of vanishing of the functions Q r (r] a 2 ,y) 
and Z r (r;a 2 ,y) is given in [24j . [25) . Here we present only an updated summary of 
the results. However, the radii of circular orbits where Z r (r; a 2 ,y) vanishes are closely 
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Figure 2. Classification of the KdS spacetimes according to the number of 
circular orbits where Q r (r; a 2 ,y) = (first digit) and Z r (r;a 2 ,y) = (second 
digit). In the parameter plane (a 2 ,y), the functions max (a 2 ) (upper dashed- 
dotted) and yh,min( a2 ) (lower dashed-dotted) determining the extrema of the 
function y^(r;a 2 ) separate regions corresponding to the black-hole spacetimes 
BH-22 and to the naked-singularity spacetimes. The left solid and right solid 
curves determine the naked-singularity spacetimes corresponding to the classes 
NS-32 and NS-21 and separate the naked-singularity spacetimes into regions 
corresponding to the classes NS-11, NS-21, NS-31, NS-32 and NS-33. 



related to embedding diagrams of the equatorial plane of the optical geometry and 
will be discussed in more details in section [6l 

There are one class of the KdS black- hole spacetimes and five classes of the naked- 
singularity spacetimes differing in the number of circular orbits where G r (r;a 2 ,y) 
vanishes and orbits where Z r (r;a 2 ,y) vanishes (see figure [2]). Examples of different 
types of behaviour of the functions are illustrated in figure [3] 

5. Relativistic dynamics in terms of inertial forces 

In few examples we show how the results obtained by relatively complex standard 
general relativistic methods could be easily and straightforwardly obtained in the 
framework of the inertial forces approach. Therefore we believe that the method of 
inertial forces, combined with standard approaches, could bring in an effective way 
some new results in modelling accretion disks. 

5.1. Equatorial circular motion 

Acceleration necessary to keep a particle in an uniform circular motion in the 
equatorial plane (0 — it/2) can be expressed due to relation f8} in the simple form [5], 



where — oo < v < oo, while — 1 < v < 1. This equation enables an effective discussion 
of the properties of both accelerated and geodesic motions. Remind that there are 
only radial components of the inertial forces non-vanishing and = = in the 
equatorial plane. 



a r (r, v) 



g r -i 2 z r -v{i + v 2 ) 1/2 c 



(47) 
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Figure 3. Gravitational acceleration Q r (solid) and velocity independent 
parts of the centrifugal Z r (dashed) and Coriolis C r (dotted) accelerations of 
particles moving along equatorial circular orbits in the KdS spacetimes. There 
are examples of qualitatively different types of behaviour of Q r (r;a?,y) and 
Z r (r; a 2 ,y) corresponding to the classification of the spacetimes according to the 
number of circular orbits where Q r (r;a 2 ,y) = and Z r (r; a 2 ,y) = 0. Positive 
(negative) parts of Q r (r\a 2 ,y) and Z r (r; a 2 ,y) determine repulsive (attractive) 
forces oriented outwards (towards) the singularity at r = 0. The positive function 
C r (r; a 2 ,y) determines the repulsive (attractive) Coriolis force in the case of the 
corotating orbits with v > (countcrrotating orbits with v < 0). The zero 
points of the functions determine the radii of circular orbits where the forces 
vanish independently of the velocity. The vertical lines denote the radii of event 
horizons. 
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For ultrarelativistic particles (v 2 ^> 1) we obtain from equation (|37j) the 
asymptotic equation 

a^{r,i) K-Q r T\Cr-v 2 {Z r ±C r ), (48) 

for photons (v 2 — > oo) we obtain the circular geodesic motion determined by the 
equation 

Z r ± C r = 0, (49) 

where the upper signs correspond to the corotating motion [v > 0) and the lower signs 
to the counterrotating motion (v < 0). Solving the equivalent equation Z 2 — C 2 = 0, 
we directly obtain the relation for the photon circular orbits in the KdS spacetimes 

EH 

y 2 aV + 2ya 2 r 2 (r + 3) + r{r - 3) 2 - 4a 2 . (50) 

This equation will be discussed carefully in section [6] in a relation with embedding 
diagrams of the optical geometry. The ultrarelativistic particles moving along the 
circular orbit at the radius of the photon circular orbits are kept by the acceleration 

aHr) = -Gr T \c r , (51) 

which is velocity independent with an accuracy 0(v~ 2 ). 

In static spacetimes, there is C r = 0, and at the radius of the photon circular 
orbit the acceleration of particles is 

4{r) - -g r , (52) 

and it is velocity independent exactly. Relation (|52|) is not limited to the case of 
ultrarelativistic particles, because Z r = at the radius of the photon circular orbit in 
static spacetimes. 

5. 2. Limitation of circular geodesies 

Velocities of particles moving along equatorial circular geodesies arc due to relations 
l[47|) and a^r (r) = determined by the expression 



^ \C 2 . - Q r Z r t hCr{C 2 - AQ r Z r + 4g 2 ) 1 / 2 il/2 

v ± = ±1 z 2 -c 2 J ' ^ ' 



It is clear that for the geodesical circular motion at a fixed radius, the condition 

C 2 - AG r Z r + 4G 2 > (54) 

must be satisfied. In the equatorial plane of the KdS spacetimes, this condition directly 
implies the relation 

^ > 0, (55) 

giving the limitation for the circular geodesies in the form r < y -1 ^ 3 = r s , i.e., we 
immediately obtain the limit given by the so-called static radius [22j . 

In static spacetimes, where C r = 0, the reality condition for the geodesies 
is Q r Z r < 0, i.e., the gravitational and centrifugal forces must point in opposite 
directions. 
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5.3. Static equilibrium positions 

A test particle is thought to be in a static equilibrium position, if it is at the 'rest' 
relative to the reference coordinate system, i.e., <f> = const, r = const and 9 — const, 
whereas its worldline is a geodesic, which means a^r = aj- = = 0. It is equivalent 
to the circular geodesical motion at a fixed radius and latitude, with Q = 0, i.e., 

v = v s = -RQlnrf- (56) 

Defining v s = "fv s , we obtain two equilibrium conditions 

g r + i 2 Z r + v s (l + v s 2 )^ 2 C r = 0, (57) 

Qg + ~V 2 s Zg + V s (l + V s 2 ) 1/2 Cg = (58) 

that have to be satisfied simultaneously. 

In the equatorial plane (6 = tt/2), condition (|55|) is automatically satisfied 
because of the vanishing of latitudinal components of all the forces. By using relations 
(|39p - ((lT|) . condition (|57[) implies the equation 

^=°< <»> 
i.e., the static equilibrium positions in the equatorial plane are really located at the 
static radius r s = y^ 1 ^ 3 . On the symmetry axis (9 = 0), condition (|55|) is also 
automatically satisfied and condition (|5T[) implies the equation 

a 2 - r 2 + yria 2 + r 2 ) 2 = Q 

implicitly determining the static radii on the axis of symmetry derived in a standard 
way and discussed in [55] . 

5.4- Perfect fluid toroidal configurations 

By projecting the conservation laws VfeT 4fc = 0, where T %k = (p + e)u l u + pg lk is the 
stress-energy tensor for the perfect fluid, into the 3-dimensional space of the special 
observers ((2|) by using the projection tensor hj — 5j+n l nj, and by using the condition 

Ul V k T lk = 0, (61) 

we arrive at the equation 

7 2 «r 7 (n fc + VT k )V k p + WjP + n 3 n k V k p + (p + e)F^ = 0. (62) 

Assuming now a circular motion of the baryotropic fluid, i.e., body with the equation 
of state p — p(e) and p ^ p(t, 4>), in the axially symmetric spacetimes, i.e., the vector t % 
and velocity v given by relations (|17p and (|18p . and ri 1 given by relation (|13p . we arrive 
at the relativistic Euler equation in terms of the inertial forces in the 'Newtonian' form 

V jP +(p + e)F; X = 0. (63) 

Of course, here, there are the only radial and latitudinal components of the inertial 
force F'j © non- vanishing. Then, in the case of the uniform distribution of the 
specific angular momentum I = L/E = const, the integration of equation J63J directly 
implies the equipressure surfaces, determined by the relation 

p(r 9) = In [ p 2 A r A e sm 2 9 1/2 

derived in a standard way in [20] 
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6. Embedding diagrams 

Properties of the optical reference geometry can be represented by the embedding 
of the equatorial plane into the 3-dimensional Euclidean space with the line element 
expressed in the cylindrical coordinates (p, z,a) in the form 

da 2 = dp 2 + p 2 da 2 + dz 2 . (65) 

It should be stressed that such embedding diagrams are very useful in a representation 
of some properties of the spacetime relevant in astrophysically important situations; 
especially behaviour of the centrifugal force, directly demonstrated by the embedding 
of the optical reference geometry, is crucial for understanding the equilibrium 
configurations of perfect fluid important in modelling accretion discs. Here, we shall 
follow the method developed in the case of Kerr-Newman geometry in [14j . 



6.1. Embedding formula 

The embedding diagram is characterized by the embedding formula z = z(p) 
determining a surface in the Euclidean space with the line element 

dz 2 

dl 2 {E) = ll + (-)]dp 2 + p 2 da 2 , (66) 

isometric to the 2-dimensional equatorial plane of the optical space determined by the 
line element 

dl = h rr dr + h^dcj) 2 , (67) 

where the metric coefficients of the optical geometry in the KdS background are given 
by the relations 

r(l + a 2 y) 2 [r 3 + a 4 ry + a 2 (2 + r + r 3 y)] 

r 

[r 3 + a A ry + a 2 (2 + r + r 3 y)] 2 



~ i y± -r u- yj L' T" u, i y -r u, ^ -r i ~r i yjj , , 

h rr — — , (D8J 



2 A r 



(69) 



The azimuthal coordinates can be identified (0 = a) and we can put p 2 = h^. Thus 
the differential form of the embedding formula is governed by the relation 

<!> 2 = M|)*-1 (70) 

It is convenient to transfer it into the parametric form z(p) = z(r(p)) with r being 
the parameter, i.e., 



whereas the sign in this formula is irrelevant, leading to isometric surfaces |17j . 

Using the relation p 2 — and the metric coefficients (|68|) and (|69|. we obtain 
the differential form of the embedding formula 



S^y^F' (72) 

where 

L = Ay 4 + By 3 + Cy 2 + Dy + E (73) 
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and 

A = a 6 r 6 (a 2 + r 2 ) 3 , (74) 
B = aV(a 2 + r 2 ) 2 [3r 3 + a 2 (3r + 10)], (75) 

C = -a 2 r 3 (a 2 + r 2 )[-3r 7 - a 2 (6r 2 + 16r + 9)r 3 - (76) 
a 4 r(3r 2 + 12r + 37) + 4a 6 ], 

D = r 2 {r 10 + a 2 r 6 {3r 2 + 2r - 18) + a 4 r 3 [r(3r 2 - 4r - 16) + 36] -£77) 
a 6 r[r(r 2 - Ur - 18) + 60] - 4a 8 (2r + 5)}, 

E = 4a 8 - 4a 6 r[r(r - 3) + 6] - 3a 4 (r - 2)r 2 [r(4r - 3) + 6] - (78) 
2a 2 r 5 [r(6r - 17) + 18] + (9 - 4r)r 8 . 

6. 2. Features of embedding diagrams 

6.2.1. Limits of embeddability From the differential form of the embedding formula 
(|72p. it is clear that due to the condition (dz/dr) 2 > 0, the equatorial plane of the 
optical geometry is not entirely embeddable into the 3-dimensional Euclidean space. 
The embeddable regions are determined by the condition L(r;a 2 ,y) < 0, whereas 
the equality in this relation determines the limits of embeddability. The equation 
L(r;a 2 ,y) — implicitly defines two functions yL±{r;a 2 ), whereas yL-{r;a 2 ) < 
j/L + (r; a 2 ) for all values of a 2 and r. Instead of giving long explicit expressions for the 
functions yL±{f]a 2 ), wc only present their properties and examples of different types 
of their behaviour (see figures [4] and [5J . For given values of a 2 and y, the limits of 
embeddability are determined by solutions of the equations 

y = VL±(r;a 2 ). (79) 

Because of the reality condition A r (r; a 2 ,y) > and the repulsive cosmological 
constant, the solutions are restricted by the condition 

0<y<y h (r;a 2 ). (80) 

Thus the embeddable regions are given by the condition < yL-(r;a 2 ) < y and 
V < VL+(r;a 2 ). 

The zero points of the functions yL±{f',a 2 ) are determined by the condition 
E(r;a 2 ) — (see equation (|73|) ), i.e., by solutions the equation 

4a 8 - 4a 6 r[r(r - 3) + 6] - 3a 4 (r - 2)r 2 [r(4r - 3) + 6] - 

2a 2 r 5 [r(6r - 17) + 18] + (9 - 4r)r 8 = 0, (81) 

implicitly defining functions a 2 0± (r). The function a 2 L0+ (r) has two local maxima 
a Lo+,maxi = 1-1354 and a 2 L Q + rnax2 = 1.0754 located at r = 0.6895 and at r — 1.3317. 
There is also one local minimum a 2 L0+ min = 1 located at r = 1. On the other hand, 
the function a 2 _(r) has no local extrema. 

The number of solutions of equation (|79p depends on the number of extrema of the 
functions y — yL±{i",a 2 ). They are governed by the functions a 2 Lel (r) and o? Le2 ±(r), 
which are defined implicitly by eliminating y from the equations L(r;a 2 ,y) = 
and d r L(r; a 2 , y) = 0. Since these expressions are too large to be presented here, 
we only show behaviour of the functions a 2 Lel (r) and a 2 Le2± (r) (see figure H]). The 
function a 2 Lel {r) is identical with the function a\ e (r) (see relation (|3"T)l ) governing 
the extrema of the function yh(r;a 2 ). Thus two extrema of J/L±(^;a 2 ) coalesce 
with two extrema of yh(r;a 2 ). These two extrema are the only common points of 
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?/L±(r;a 2 ) and y/i(r;a 2 ), and there is 2/L±(? , ;a 2 ) < y/j(r;a 2 ) except these points. The 
maximum of the function a 2 Lel (r) takes the value a 2 Lel rnax = 1.2120 and is located 
at r = 1.6160. The common point of a 2 Lel (r) and a 2 L0+ (r) coalesces with the local 
minimum of a 2 0+ (r) and divides the function a? Lel {r) into two parts governing positive 
(relevant for the classification) and negative (irrelevant) extrema of yz,±(r; a 2 ). The 
function a 2 Le2+ {r) has local maximum a 2 Le2+ max = 1.4706 located at r = 1.0961 and 
local minimum a 2 Le2+ min = 1.0683 located at r = 1.3289. Two of three common 
points of a 2 Le2+ {r) an d o 2 La+ {r) coalesce with local maxima of a 2 L0+ (r) and divide the 
function a 2 Le2+ (r) into parts governing positive and negative extrema of yL±{r\a 2 ). 
The function a 2 Le2 _(r) completely governs negative local extrema of yL±{r; a 2 ), which 
we do not consider here. 

All the characteristic functions a 2 j0± (r), a| el± (r) and a 2 Le2 (r) are illustrated in 
figure 2] These functions enable us to understand the behaviour of the functions 
VL±{r;a 2 ) and classify the KdS spacetimes according to the number of embeddable 
regions. 



6.2.2. Turning points of embedding diagrams In general, the plasticity of the 
embedding diagrams is characterized by the different steepness and by the turning 
points of the diagrams. Clearly, the turning points of the embedding diagrams, i.e., 
their bellies and throats, are determined by the relation dp/dz = 0. Because of the 
relation 

dz = dzdr ^ 
dp dr dp'' 

we have the condition dp/dr — 0. Due to relations ([4|, (|23|) and p 2 — h^, we can 
write dp/dr = dR/dr = 0. Considering relation ([25]) . it is clear that the turning points 
of the embedding diagrams are located at the radii where the velocity independent 
part of the centrifugal acceleration vanishes. 
Since there is 

J = {r 2 Ap}- VV(a 2 + r 2 )y 2 + r 2 a 2 [(2r + 5)a 2 + r 2 (2r + 3)}y + 
r 4 (r-3) + ra 2 [r(r-3) + 6] -2a 4 }, (83) 

then for given values a 2 and y, the radii of turning points of the embedding diagrams 
are given by solutions of the equation 

y = y T (r; a 2 ) = {2a 2 r 2 (a 2 + r 2 )}" 1 {-(2r + 5)ra 2 - r 3 (2r + 3) + 

yjr{a 2 + 3r 2 )[8a 4 + m 2 (16r + 1) + r 3 (8r + 3)]} (84) 

restricted by the condition 

< y L -(r;a 2 ) < y < y L+ (r;a 2 ). (85) 

The zero points of the function yT(r;a 2 ) are determined by solutions of the 
equation 

r 4 (r - 3) + ra 2 [r(r - 3) + 6] - 2a 4 = 0, (86) 

which we consider as an implicit form of the functions ay 0± (r). The maximum of this 
function takes the value a^ 0± max = 1.3667 located at r = 0.8116. 
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Figure 4. Characteristic functions af,o±( r ) (closely dotted) governing zeros of 
the function yL±{ r \ a2 )\ a z, e i( r ) (dashed-dotted) and a\ e 2±{ r ) (solid) governing 
extrema of J/z,±(r; a 2 ); a™i (r) (dotted) governing zeros of the function i/j^rja 2 ); 
a Te±( r ) (dashed) governing extrema of i/t(x; a 2 ). The thick parts of the dashed- 
dotted and solid curves govern positive extrema of yL±(r \ a 2 ). In the parameter 
line (a 2 ) of the spacetimes, the extrema and common points of the characteristic 
functions separate regions corresponding to different types of the behaviour of 
DL±(r;a 2 ) (see figure 15) , 



The number of solutions of equation (|84[) depends on the number of the extrema 
of the function yT{r\a 2 ). They are determined by solutions of the equation 

{^a 4 r 3 (a 2 + 3r 2 )[8a 4 + r(l6r + l)a 2 + r 3 (8r + 3)]}" 1 x 

{r[12a 8 + r(48r + l)a 6 + 3r 3 (24r + l)a 4 + 3r 5 (16r + l)a 2 + 

3r 7 (4r + 3)]a 2 } - 5a 4 - 12rV - 3r 4 = 0, (87) 

which we consider as an implicit form of the functions a^ e _(_(r) with the maximum 
a Te+,max = 1-8126 located at r = 1.3285. 

The characteristic functions a^j- (r) and a™± ( r ) are illustrated in figure 3] and 
enable us to understand the behaviour of the functions yr(r; a 2 ) and classify the KdS 
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Table 1. Different types of the KdS spacetimes determined by corresponding 
different types of behaviour of the functions J/L±(r;a 2 ) and j/r(r;a 2 ). The 
number of extrema y L+iel (a 2 ), VL+,e2{a 2 ), yL+,e3(a 2 ), VL+,e4.{a 2 ), yL~,max{a 2 ), 
yT,min{a 2 ) and yT,max(a 2 ) of the functions yL±{r;a 2 ) and y T (r\a 2 ) satisfying 
the condition < y < yh(r;a 2 ) are subsequently expressed by the digits in the 
last column. The types Ea, Eb and Ec differ in the relation among the values of 
extrema yL+,e2( a2 )y 2/L+,e3( a2 ) an d yh+^iifl 2 ) an d the types Fa and Fb differ 
in the relation among the values of extrema yL+^^a 2 ) and yT,max(a 2 )- 



T 


Range of a 2 


Extrema 


T 


Range of a 2 


Extrema 


A 


(0,1) 


0,1,0,0,0,0,0 


Ec 


(1.2032, 1.2120) 


1,1,1,1,0,0,0 


SO 


1 


0,1,0,0,0,0,0 


S2 


1.2120 


0,0,1,1,0,0,1 


B 


(1, 1.0683) 


1,1,0,0,0,0,0 


Fa 


(1.2120, 1.2938) 


0,0,1,1,0,0,1 


C 


(1.0683, 1.0754) 


1,1,0,1,1,0,0 


Fb 


(1.2938, 1.3667) 


0,0,1,1,0,0,1 


D 


(1.0754, 1.1354) 


1,1,0,1,0,0,0 


G 


(1.3667, 1.4706) 


0,0,1,1,0,1,1 


Ea 


(1.1354, 1.1597) 


1,1,1,1,0,0,0 


H 


(1.4706, 1.8126) 


0,0,0,0,0,1,1 


SI 


1.1597 


0,1,1,0,0,0,0 


1 


(1.8126, oo) 


0,0,0,0,0,0,0 


Eb 


(1.1597, 1.2032) 


1,1,1,1,0,0,0 









spacetimes according to the number of turning points of the embedding diagrams. 

6.2.3. Classification The number of embeddable regions and turning points of the 
embedding diagrams is determined by the number of solutions of equations (|79p 
and (|84j) . which depends on the number of extrema of the functions yL±(r;a 2 ) and 
yT{r;a 2 ). Therefore we denote yL-,max(,a 2 ) as the only maximum of y£_(r;a 2 ) 
which becomes positive (relevant for the classification) for some values of a 2 and 
f, VL+.eiio 2 ), UL+,e2(a 2 ) as the extrema of yL+(r;a 2 ) coalescent with the minimum 
and maximum of the function yh{r;a 2 ); yL+,es{o, 2 ), yL+,e4(a 2 ) as the remaining two 
extrema of y£+(r; a 2 ). Finally we denote yT,min(a 2 ) and yT,max{a> 2 ) as the minimum 
and maximum of the function y-rf^; a 2 )- 

Using the characteristic functions a^ e± (r), a^ 0± (r), a 2 L0± (r), a 2 Lel {r) and 
a ie2±( r )i we can distinguish different types of behaviour of the functions yL±(r;a 2 ) 
and yT(r;a 2 ) differing in the number of their extrema satisfying the condition 
< y < yh(r;a 2 ) (see figure [5]). The different types of the KdS spacetimes obtained 
this way are shown in table [T] The classification of the KdS spacetimes according to 
the number of embeddable regions and turning points of the embedding diagrams can 
be now given in the following way. We step by step discuss the number of embeddable 
regions using the number of solutions of equation (f79|) restricted by condition (|80p for 
given values of a 2 and y, and the number of turning points of embedding diagrams 
given by the solutions of equation (f84|) . restricted by condition (|85]) . for each type of 
the KdS spacetimes separately. 

In the spacetimes of type 'A', for < y < yi+, e 2(a 2 ), there are four solutions 
of equation (|79|) satisfying condition (|80|) and then two embeddable regions. There 
are two solutions of equation (|84j) satisfying condition ()85[) . i.e., two turning points 
of the embedding diagrams. Moreover, there are three solutions of equation ()35|1 
determining three even horizons and thus two stationary regions of the spacetimes. 
We denote the related class as BH22, alias the KdS black-hole spacetimes with two 
embeddable regions and two turning points. For y > y_L+,e2(a 2 )j there are two 
solutions of equation (I79| satisfying condition (|80|) and thus only one embeddable 
region. The only solution of equation (|84|) satisfying condition (|85|) determines one 




Figure 5. Functions yL±( r ', a2 ) (solid) limiting cmbeddablc regions (white); 
2/-r(r;a 2 ) (dashed) determining turning points of embedding diagrams; j/^(r;a 2 ) 
(dashed-dotted) determining locations of event horizons and limiting the dynamic 
regions (gray). For illustration, the behaviour of the function y p /i(r;a 2 ) (dotted) 
governing radii of photon circular orbits is also given. For given a 2 and y > 0, 
limits of embeddability are determined by solutions of y = yL±(r; a 2 ) that satisfy 
the condition y < j/^(r;a ). Turning points of the embedding diagrams are 
determined by solutions of y = yj-{r\a ) and cmbeddablc photon circular orbits by 
solutions of y = y ph (r; a 2 ), whereas the solutions must satisfy the embeddability 
condition j/£_(r; a 2 ) < y < yL+(r; a 2 ). In the parameter line (y) of the KdS 
spacctimes, extrema of 2AL±(r;a 2 ) and yj-(r\a 2 ) separate regions corresponding 
to different classes of the spacetimes differing in the number of cmbeddablc regions 
(first digit) and turning points of the embedding diagrams (second digit). Note 
that in the case of the classes NS12 and NS22, the second digits exceptionally 
denote one turning point and one inflexion point of the diagrams. 
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Figure 5. (Continued.) 
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Figure 5. (Continued.) 
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Table 2. Classification of the KdS spacetimes according to the number of 
embeddable regions (first digit) and turning points of the embedding diagrams 
(second digit). Limits of the ranges of the parameter y (being functions of the 
parameter a 2 ) are illustrated in figure [8] 



Class 


Range of y 


Class 


Range of y 


Type A, 


SO 


Type Ec -continued 


BH22 


(0;yL+M a ^)) 


BH32 


(yL+,ei(a. );yL+,e2(a )) 


NS11 


(VL+,e2(a 2 ); OO) 


NS21 


{yL+,e2(a );y i+ , e 4(a )) 


Type B 




IN oil 


(yL+,ei(a ); oo) 


NS33 


(0;yL+, e iK)) 


lype oz 




BH22 


(yi + ,el(ct 2 ); VL+,e2(a 2 )) 


NS13 


(ci ( 'AW 

(0;yL+,e3(a )) 


NS11 
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NS23 




Type C 
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{yL + ,el(a. 2 );yL+,e4(a 2 )) 


NS13 


(yT,min{a 2 );yT,max(a 2 )) 


BH22 


(yL+,ei(a. 2 );y L+>e2 (a 2 )) 


NS12 


y = yT,max(a 2 ) 


NS11 


(yL+,e2(a. 2 ); OO) 


NS11 


{VT,max(a 2 ); oo) 


Type Ec 




Type I 




NS13 


(0;yL+M a2 )) 


NS11 


(0;oo) 


NS23 


{yL+,e3{a 2 );yL+,ei(a 2 )) 







turning point of the embedding diagrams. Moreover, there is one solution of equation 
(135 j) . i.e., one cosmological horizon separating stationary region and dynamic region. 
We denote the related class as NS11, alias the KdS naked-singularity spacetimes with 
one embeddable region and one turning point. Note that although there are two 
horizons and two stationary regions in the case of y = z/L+,e2(ft 2 ), the outer stationary 
region coincides with the outer horizon where the geometry is not embeddable and 
thus we cannot consider it. Using the same way, we can sort the remaining types of 
the spacetimes obtaining the classification for all the KdS spacetimes (see table [2]) . 

Now it is clear that the functions y L +, e i{a 2 ), VL+.e2(a 2 ), VL+,eii.a 2 ), 2/L+,e4(a 2 ), 
VL-, m ax{a 2 ), VTjnin{a 2 ) and yT,max{a 2 ) separate the parametric plane {a 2 ,y) into 
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regions corresponding to different classes of the KdS black-hole spacetimes and naked- 
singularity spacetimes differing in the number of embeddable regions and number of 
turning points of the embedding diagrams (see figure [S]). Qualitatively different types 
of the embedding diagrams corresponding to the presented classification are illustrated 
in figures [9]fT0l 

6.3. Embeddability of photon circular orbits 

The radii of photon circular orbits are given by relation (|49| . Thus, of course, in the 
stationary (non-static) spacetimes, where C r ^ 0, the radii of photon circular orbits 
do not coalesce with the radii of the turning points of embedding diagrams (radii of 
circular orbits where Z r — 0), which holds for the static spacetimes, where C r = 0. 

For given values of a 2 and y, radii of circular photon orbits are determined (due 
to equation ([50)0 by solutions of the equation 

y = y ph (r;a 2 ) = * '- (88) 

restricted by condition (|36[) . The zero points of the function are determined by 
solutions of the equation 

a 2 = a% (r) = ±(r-3)\. (89) 

The number of circular orbits is determined by the extrema of the function, which 
coalesce with the extrema of the function j/^(r;a 2 ). Thus they are governed by the 
function af le (r) given by equation ([57)) . Solutions of equation ([55)1 satisfying condition 
(|85[) then correspond to the embeddable orbits, number of which depends on the 
number of common points of y p h{r; & 2 ) and yL±{f", a 2 ). The points are determined by 
solutions of the equation 

[a 4 + (r - 3)r 3 + a 2 r(2r - 1 )]{16a 10 + 16a 8 (23 - 4r)r 2 - 

9r 8 (4r - 3) 3 + 4aV[9 + 32(45 - llr)r] - 8aV(8r(8r - 37) - 9] + 

3a 2 r 6 [27 - 16r[9 + r(32r - 105)]]} = 0, (90) 

implicitly defining the function a 2 hL (r), which determines a first kind of common 
points, and the function a 2 e (r), determining the other kind, coalescing with the 
extrema of yh(r;a 2 ). All the characteristic functions a 2 hL (r), af le (r) and a 2 h0 (r) are 
illustrated in figure El 

6.3.1. Classification Number of photon circular orbits depends on the number of 
extrema of the function y p h(r; a 2 ). Since these extrema coalesce with the extrema of 
the function yh{r; a 2 ) and there is y p h{r; a 2 ) < yh(r; a 2 ) for all the values of r > and 
a 2 , it clear that the KdS black-hole spacetimes contain always two photon circular 
orbits in the outer stationary region and one orbit in the inner one. In the KdS 
naked-singularity spacetimes, only one such orbit always exists. However, all the 
orbits do not occur in the embeddable regions. The number of embeddable photon 
circular orbits depends on the number of common points of the functions y p h(r; a 2 ) and 
y £± (r;a 2 ) and their relative positions. We denote y p h,L(a 2 ) as the value of the function 
y p h{r; a 2 ) corresponding to the intersection of the functions y p h(r] a 2 ) and yL±{r; a ), 
determined by the function a 2 hL (r), i.e., those intersections, which do not correspond 
to the extrema of the function y p h(r;a 2 ). By using the characteristic functions 
a phL ( r ) > a he ( r ) an< l a pho ( r ) ' we can distinguish different types of the behaviour of the 
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Figure 6. Characteristic functions a? hIj (r) (solid) governing common points 
of the functions y ph (r;a 2 ) and ?/L±(»*;a 2 ); (r) (dashed-dotted) governing the 
extrema of y p h(i", a 2 ) coalescent with the extrema of the function y^(r;a 2 ) and 
with the other common points of y p h{r;a 2 ) and j/L±(r; a 2 ); a ph0 (r) (dotted) 
governing zeros of y p h(r; a?). The physical relevant solutions are then denoted by 
the thick parts of the relevant curves. In the parameter line (a 2 ) of the spacetimes, 
the common points of the characteristic functions separate regions corresponding 
to the different types of behaviour of y p h{r; a2 ) ( see figure [7J. 



Table 3. Different types of the KdS spacetimes determined by the corresponding 
different types of behaviour of the functions yL±{ r \ a ' 2 ) an d t/ p h(r;a 2 ). The 
number of common points (except the common points coalescent with the extrema 
of 2/?,(r;a 2 )) is expressed by the digits in related column, whereas the first digit 
denotes the number of common points 'between' the extrema of y p h{r; ° 2 ) an d the 
second one denotes the number of extrema 'behind' the maximum of y ph (r; a 2 ). 
Note that the types si and s2 represent limit cases of the behaviour of the 
functions. There is one negative common point of the functions in the case of 
the type 'a', denoted by the negative sign at the related digit. 



T 


Range of a 


Common points 


a 


(0,0.5057) 


-1,0 


si 


0.5057 


0.0 


b 


(0.5057, 0.9718) 


1,0 


s2 


0.9068 


0.0 


c 


(0.9068, oo) 


0.1 



functions y p h(r;a 2 ) and yL±{r' 1 a 2 ) (see figure [7] and table [3|). In the parameter line 
(y) of the KdS spacetimes, the common points of y p h(r;a 2 ) and yL±{r;a 2 ) separate 
regions corresponding to different classes of the KdS spacetimes with two, one ore 
none embeddable photon circular orbits (see figure [7]) . The classification can be then 
obtained in the same way as the previous ones, i.e., by step by step discussing the 
number of embeddable circular photon orbits for each type of the spacetimes separately 
(see table |4]). The different classes of the KdS spacetimes are illustrated in figured! 





Figure 7. Functions ?/ph(r;a 2 ) (dotted) governing the radii of photon circular 
orbits; 2AL±(r;a 2 ) (solid) limiting embeddable regions (white); yj<(r;a 2 ) (dashed) 
determining turning points of embedding diagrams; Jfft(r;a 2 ) (dashed-dotted) 
determining locations of event horizons and limiting the dynamic regions (gray). 
In the parameter line (y) of the KdS spacetimes, the common points of 
Vphii", ?) an d VL±{'>",a 2 ) separate regions corresponding to different classes of 
the spacetimes differing in the number of embeddable photon circular orbits. 



Table 4. Classification of the KdS spacetimes according to the number of 
embeddable photon circular orbits. Limits of the ranges of the parameter y (being 
functions of the parameter a 2 ) are illustrated in figure [8] 



Class 


Range of y 


Class 


Range of y 


Type a 


si 


Type s2 




-2 




-1 


(0;yh,m,ax(a 2 )) 


-0 


(yh,max(a. 2 ); oo) 


-0 


(ft,m«i(« 2 ); OO) 


Type b 




Type c 




-1 


(0;y phi (a 2 )) 


-1 


(0;y phL (a 2 )) 


-2 


(yphL(.a 2 );yh, m ax(a 2 )) 


-0 


(y P hL(a 2 ); oo) 


-0 


(yh,max(a 2 ); oo) 
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Figure 8. Classification of the KdS spacetimes according to the number of 
cmbcddable regions (first digit), turning points of embedding diagrams (second 
digit) and embeddable photon circular orbits (digit following dash). In the 
parameter plane (a 2 ,y), the functions yL+,el( a2 ) an d UL+,e2( a2 ) (dashed- 
dotted) governing the extrema of the function yL±(r;a 2 ) coalescent with 
extrema of the function y h (r;a 2 ) (and y ph (r; a 2 )); yL+, e :i( a2 ), yL+,e4.(a 2 ) and 
yL-,max{a 2 ) (solid) governing the remaining extrema of yL±(r; a 2 ); y p h,L(a 2 ) 
(dotted) governing the remaining common points of the functions y p h( r ', a2 ) 
and i/£±(r;a 2 ); yT,min( a2 ) an d tlT,max( a2 ) (dashed) governing extrema of the 
function yT(r;a 2 ) separate regions corresponding to the different classes of the 
KdS spacetimes. Note that in the case of the classes NS12-1 and NS22-0, the 
second digits exceptionally denote one turning point and one inflexion point. 
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Figure 9. Embedding diagrams of the optical reference geometry of the KdS 
spacctime of the classes BH22 — 1 and NS23 — 1. 



7. Conclusions 

The presented study of the optical reference geometry and related inertial forces in the 
Kerr-de Sitter (KdS) spacetimes represents a generalization of results obtained earlier 
for the simpler limit cases of the KdS spacetimes, i.e., for the Kerr and Schwarzschild- 
de Sitter (SdS) spacetimes (see [T3], [5] and [S])- 

In the KdS spacetimes, effects of the rotation of source and the cosmic 
repulsion combine in a relatively complex way, being described by the rotational and 
cosmological parameters a 2 and y. In order to better understand the interplay of both 
the parameters, it is useful to summarize results of our investigation and compare 
them just with the studies on the background of the Kerr, SdS and Schwarzschild 
spacetimes, where the rotational and cosmological effects manifest separately. 

The KdS black-hole spacetimes admit two stationary regions and two dynamic 
regions, separated by two black-hole horizons and one cosmological horizon. The KdS 
naked-singularity spacetimes then contain one stationary and one dynamic region 
separated by the cosmological horizon. Note that the optical reference geometry 
and related inertial forces are defined in the stationary (static) regions only. The 
stationary and axially symmetric Kerr black- hole spacetimes (a 2 < 1, y = 0) admit 
two stationary regions and one dynamic region, separated by two black-hole horizons. 
In the Kerr naked-singularity spacetimes (a 2 > 1, y = 0), only one stationary region 
remains. 
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Figure 10. Embedding diagrams of the optical reference geometry of the KdS 
spacctime of the classes NS13 — 1 and NS12 — 1. 

The static and spherically symmetric SdS black-hole spacetimes (a 2 = 0, < 
y < 1/27) contain one static region, and two dynamic regions, separated by black-hole 
and cosmological horizons. In the SdS naked-singularity spacetimes (a 2 = 0, y > 
1/27), only one dynamic region exists. Therefore we do not pay attention to these 
naked-singularity spacetimes. The well known Schwarzschild black-hole spacetimes 
(a 2 = 0, y = 0) contain the only horizon separating one dynamic region and one 
static region. 

An investigation of the optical geometry in a given spacetime can be considered 
divided into parts concerning the behaviour of the inertial forces, visualization of the 
optical geometry itself via embedding diagrams, and application of the formalism in 
relativistic dynamics. 

Of course, compared to the simpler cases of the Schwarzschild, SdS and Kerr 
spacetimes, the behaviour of the forces is more complex in the KdS spacetimes, 
incorporating the influence of both the spacetime parameters on the particle motion. 
There are more points of their vanishing, changing orientation and extrema. We have 
classified the KdS spacetimes according to the number of equatorial circular orbits 
where the velocity independent gravitational force vanishes and the orbits where 
the centrifugal force vanishes independently of the velocity (see table [SJ . We can 
summarize that 

• KdS black-hole spacetimes contain only one circular orbit where the gravitational 
force vanishes, and only one orbit where the centrifugal force vanishes 
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Table 5. Classification of the KdS, Kerr, SdS and Schwarzschild spacetimes 
according to the number of circular orbits where the gravitational force vanishes 
(first digit), and orbits where the centrifugal force vanishes independently of the 
velocity (second digit). 





BH spacetimes 


NS spacetimes 


KdS 


BH-22 


NS-11, NS-21, NS-31, NS-32, NS-33 


Kerr 


BH-12 


NS-01, NS-11, NS-21, NS-22, NS-23 


SdS 


BH-11 




Schw. 


BH-01 





independently of the velocity in each of both the stationary regions. Both the 
gravitational and centrifugal forces change their orientations on these orbits. The 
same number of the orbits appears in the Kerr black-hole spacetimes, except for 
the outer stationary region, where there is no circular orbit where the gravitational 
force vanishes. In the SdS black-hole spacetimes, the gravitational force vanishes 
on the circular orbit of the static radius r — y -1 ^ 3 only. The centrifugal force 
vanishes independently of the velocity (and y) on the orbit of the radius r = 3. As 
well as in the Schwarzschild spacetimes, this is the only orbit where the centrifugal 
force vanishes independently of the velocity. On the other hand, the gravitational 
force does not vanish in the Schwarzschild spacetimes. 

• KdS naked-singularity spacetimes can contain even three circular orbits where 
the gravitational force vanishes and other three orbits where the centrifugal 
force vanishes independently of the velocity, which indicates a relatively complex 
structure of these spacetimes as a result of mixed influence of the rotation of the 
source and the cosmological repulsion. The Kerr naked-singularity spacetimes 
contain at maximum only two circular orbits where the gravitational force 
vanishes, but there are at most three circular orbits where the centrifugal force 
vanishes independently of the velocity in accord with the case of the KdS naked- 
singularity spacetimes. With the parameters a 2 and y growing, the structure of 
the KdS spacetimes simplifies. There are only one orbit with the vanishing of the 
centrifugal force and the other with the vanishing of the gravitational force. 

The velocity independent part of the Coriolis force does not vanish neither in the KdS 
nor in the Kerr spacetimes, while in the SdS and Schwarzschild spacetimes the force 
does not occur at all. 

In general, we can summarize, in static and spherically symmetric spacetimes, 
all the geodesies of the optical reference geometry coincide with trajectories of light, 
which explains the name 'optical geometry'. We have seen that test particles moving 
along the circular photon orbits are kept by the velocity-independent gravitational 
force only, because the remaining centrifugal force vanishes here independently of the 
velocity. Unfortunately, in the stationary (non-static) spacetimes, the photon circular 
orbits in the equatorial plane do not correspond to the circular trajectories where the 
centrifugal force vanishes independently of the velocity. 

As for the embedding diagrams, the equatorial plane of the optical reference 
geometry can not be entirely embeddable into the 3D Euclidean space and the 
embedding diagrams then consist from several separated parts. We have characterized 
the shape of the embedding diagrams by the number of their turning points coalescing 
with the radii of circular orbits where the centrifugal force vanishes independently 
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Table 6. Classification of the KdS, Kerr, SdS and Schwarzschild spacetimes 
according to the number of embeddable regions (first digit), turning points of 
the embedding diagrams (second digit) and the number of embeddable photon 
circular orbits (digit following dash). 





BH spacetimes 


NS spacetimes 




KdS 


BH22-0, BH32-0 


NS11-0, NS21-0, NS22-0, NS23-0, 


NS11-1, NS12-1 




BH22-1, BH22-2 


NS13-1, NS21-1, NS22-1, NS23-1, 


NS33-1 


Kerr 


BH22-2, BH22-1 


NS11-1, NS12-1, NS13-1, NS23-1, 


NS33-1 


SdS 


BH11-1 






Schw. 


BH11-1 







of the velocity. Because the radii of photon circular orbits in the equatorial are not 
located at the radii where the centrifugal force vanish, i.e., at the radii of the turning 
points of the diagrams, we have discussed embeddability of the photon circular orbits 
as well. Thus the KdS spacetimes can be also classified according to the features of 
the embedding diagrams (see table [6])). We can summarize that 

• KdS black-hole spacetimes contain three or two separated embeddable regions in 
contrast to the Kerr black-hole spacetimes, where only two regions occur, and 
to the SdS and Schwarzschild black-hole spacetimes with the only embeddable 
region. 

• KdS naked-singularity spacetimes contain one, two or three separated embeddable 
regions as well as the Kerr naked-singularity spacetimes. 

• KdS black-hole embedding diagrams contain two turning points as well as the 
Kerr black- hole diagrams. But the SdS and Schwarzschild black- hole embedding 
diagrams contain one turning point, coalescing with the radius of the circular 
orbit where the centrifugal force vanishes independently of the velocity, i.e., with 
the radius r = 3. 

• KdS naked-singularity embedding diagrams contain one, two or three turning 
points as well as the Kerr naked-singularity diagrams. 

• KdS black-hole spacetimes contain three photon circular orbits, whereas two of 
them occur in the outer stationary region and one of them in the inner one, as well 
as in the Kerr black-hole spacetimes, but in contrast to the SdS and Schwarzschild 
black-hole spacetimes, where only one photon orbit occurs at the radius r = 3. 

• KdS naked-singularity spacetimes contain one photon circular orbit as well as the 
Kerr naked-singularity spacetimes. 

• KdS black-hole spacetimes contain two, one or none embeddable photon 
circular orbits in the outer stationary region in contrast to the Kerr black-hole 
spacetimes, where always two or one embeddable orbits exist, and to the SdS 
and Schwarzschild black-hole spacetimes, where always one embeddable orbit 
occurs at the radius coalescing with the radius of turning point of the embedding 
diagram. 

• KdS naked-singularity spacetimes contain one or none embeddable photon 
circular orbit in contrast to the Kerr naked-singularity spacetimes, where always 
one embeddable orbit occurs. 
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Both the presented parts of the investigation of optical reference geometry formal- 
ism and related classification of the KdS spacetimes can be consider to be the basic 
and necessary study, before the formalism is applied in the dynamics. 

By using some particular examples of the particle and fluid orbital motion, 
we have just directly demonstrated that the optical reference geometry and related 
inertial forces formalism presents certain, we can say 'pseudo-Newtonian' approach 
in relativistic dynamics. This approach is based on the Newtonian way of dynamics 
description, but taking the general relativistic effects fully into account. Now, we can 
state that the simplicity and applicability of the formalism, which were only predicted 
in a general way in some papers (e.g. [5]) or directly shown, but only for the simpler 
types of spacetimes, 'survive' even in the complex types of the KdS spacetimes, where 
its application is then more useful all the more just because of the complexity of the 
spacetimes. 

In spite of the more or less technical aspect of the paper, the study can be reviewed 
from the astrophysical point of view as well. We have shown that the combined effect 
of rotation and cosmic repulsion is relevant in both the cases of black- hole and naked- 
singularity backgrounds. Clearly, choosing a characteristic type of motion and knowing 
the related behaviour of the forces, we can understand the influence of the rotational 
and cosmic parameters. Of course, there are many astrophysical and conceptual 
reasons to choose the circular motion, as we have done. Among others, comparing the 
results for the Kerr and KdS spacetimes, the additional radius in the KdS spacetimes, 
where the gravitational force vanishes, the same number of circular orbits where the 
centrifugal force vanishes and qualitatively same types of behaviour of embedding 
diagrams suggest that the influence of the cosmological constant is fundamentally 
imprinted mainly in the gravitational force. Note that it was predicted in [H], by 
comparing the behaviour of the gravitational force in the SdS and Schwarzschild 
spacetimes (see tables [5] and [6]) . 

One of the important new results concerns properties of the so-called static radii, 
where the gravitational attraction of a black hole (naked singularity) is just balanced 
by the cosmic repulsion. We have found that in the equatorial plane the location of 
the static radius r = y^ 1 ^ 3 is independent of the rotational parameter a, and formally 
coincides with the results from the SdS black-hole spacetimes, while on the axis of 
symmetry it depends on a. Physically, it is important that these results are easily 
implied by the analysis of the behaviour of the inertial forces. 

From the physical point of view, the most important new result is given by the 
possibility to derive the Euler equation easily in the framework of optical geometry 
and related inertial forces as shown in Section 3. It enables us to treat behaviour 
of perfect fluid orbiting black holes, and in general all the accretion processes, in the 
pseudo-Newtonian way related to the inertial forces. We believe the pseudo-Newtonian 
approach could be useful in modelling self-graviting discs around SdS and KdS black- 
holes. It should be noticed in this connection that character of the centrifugal force, 
which is relevant in determining the structure of thick discs, is clearly reflected by the 
embedding diagrams of the optical geometry (see [T4]1. 

Thus we can state, we dispose an alternative method to the standard general 
relativistic approach, practically and intuitively applicable in the KdS spacetimes, and 
hence suitable for the description of astrophysical phenomena around supermassive 
black holes in centre of giant active galactic nuclei. 
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